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Introduction 


Studies of the ecology and management of Howellia aquatilis in the Swan Valley 
have highlighted pond hydrology as a key factor in meeting the habitat needs of the 
species (Lesica, 1990; Lesica, 1992). Understanding controls on the hydrologic 
functioning of H. aquatilis ponds requires investigation of water balances and local 
hydrogeology of a representative sample of ponds. At the same time, the ecology of the 
species implies an important role for pond morphometry in the response of populations to 
changing water balances (Shapley and Lesica, 1997). Size, shape and orientation may 
also be important to other ecological factors such as light regime. Pond geometry is 
therefore viewed as one important characteristic in selecting a subsample of the more than 
100 known inhabited ponds for future investigation. The purpose of this study is to 
characterize key morphometric and geochemical variables for a large number of occupied 
ponds, addressing certain components of the ‘extensive’ data collection suggested by 
Shapley and Lesica (1997). This broad survey of habitat hydrology will provide an 
important basis for choosing sites for future site specific studies, studies needed to answer 
questions regarding source contribution and land management implications. 


Study Sites 


Howellia aquatilis was known from 101 wetland basins at the time of this survey, 
spanning the 50-mile length of the Swan River Valley in Lake and Missoula Counties, 
Montana. During initial field visits timed to coincide with active growth and 
reproduction (between July 21 and August 11, 1997) specific conductance, pH and 
temperature data were collected from all known occupied ponds on public land and from 
ponds on private land where access was permitted (n = 97). Observations on the presence 
of Howellia aquatilis and the invasive grass Phalaris arundinaceae were made incidental 
to this first 1997 visit. Pond water levels were marked for later comparison and 
determination of water level and specific conductance changes. 


A subset of the ponds was chosen for cross-sectional elevation surveys and 
morphological modeling. These ponds were selected on the basis of three criteria: 
1) Distribution among the recognized clusters of H. aquatilis occurrences within the 
Swan Valley; 
2) Gross local relief of the immediate pond watershed, based upon 1:24,000 topographic 
mapping, and 
3) Specific conductance of pond waters, measured during a comprehensive survey of 
known occupied ponds during the first 1997 site visit. 
For each recognized cluster of H. aquatilis occurrences, ponds representing high and low 
relief settings and high and low specific conductances (relative to the range displayed by 
the occurrence cluster) were selected for elevation surveys. (Relief categories were based 
upon the number of 20 foot contour intervals between the site and adjacent mapped 
drainage features, i.e., the height of local drainage divides relative to each pond.) 


Additional ponds were then chosen by lot from among the remaining occurrences in the 
cluster, up to a total approximately representative of the cluster’s contribution to the total 
set of Swan Valley occurrences of H. aquatilis. The proportionality among clusters is 
only approximate, as the time necessary for surveying individual ponds varied widely and 
was difficult to predict a priori. Thirty-five occupied sites were selected through this 
process. One site selected was an unusually large wetland that could not be effectively 
surveyed with the procedures used. Deep water prevented the completion of cross 
sections at three Condon Creek sites, resulting in under-representation of the Condon 
Creek cluster. At the remaining 31 ponds, two to four elevation cross-sections were 
completed. 


Methods 


Temperature, pH and specific conductance measurements were collected using 
one or more of three field instruments: an Extech Oyster meter with combined 
pH/conductivity electrode, an Orion Model 124 conductivity/temperature meter, and an 
Oakton Instruments pen-style pH meter. The Extech instrument was calibrated to 
multiple pH buffers on several occasions daily while in use. The performance of the 
Orion probe was checked against conductivity standards prior to each use. Under humid 
conditions the Orion unit performed erratically, and only the Extech instrument was used 
for specific conductance measurements. Agreement between the two pH meters was 
generally within 0.1 units. Where agreement was poorer, data tables present 
measurement from the calibrated (Extech) instrument. 


We measured relative elevations along surveyed cross-sections using a tripod- 
mounted surveying level and stadia rod provided by the Flathead National Forest. We 
measured linear distances with a 30 meter nylon tape. We measured elevations at 
intervals ranging from <1 meter to 10 meters, depending on pond size and pond-bottom 
relief. Stadia rod positions were usually reached on foot; on several occasions deep water 
required the use of either a kayak or an extender attached to the rod to complete cross - 
sections. The elevation of the water surface and of features indicative of full-pool 
elevation were noted along each cross - section. Compass bearings of cross sections and 
other significant landmarks are based on measurements made with a hand-held Brunton 
compass. Elevation surveys were carried out from August 27 to September 10, 1997, a 
schedule planned to coincide with seasonally low water levels and curtailed growth and 
reproduction of H. aquatilis. However, high 1997 water levels meant that many ponds 
remained far above minimum water levels during these surveys, complicating 
measurements and preventing the completion of some cross-sections. 


At each site, one cross-section followed the major axis of the pond. Additional 
sections were measured along transverse axes. For ponds with complex forms, we 
measured as many as two additional sections determined by pond geometry. We also 
attempted to map pond perimeters using a global positioning system (GPS), but found 
signal reception to be very poor at many sites due to dense forest canopy cover. 


Therefore, perimeter locations are based on the cross-section positions and occasional 
point measurements accompanied by aerial photos and field sketches of pond outlines. 
During these surveys, we examined pond perimeters for evidence of surface water inflow 
or outflow, and remeasured the pH, specific conductance and temperature of pond waters. 
The distribution of H. aquatilis and other aquatic plants along the surveyed cross-sections 
was noted qualitatively. 


Most elevations could be read with a precision of + .015 m or better, and 
elevations of the water surface are usually replicated within .01 m on a given cross 
section. However, tussocky vegetation, soft substrate and deep water substantially | 
decreased the accuracy of many measurements. The position of full-pool features 
(shorelines and vegetation transitions) is a matter of judgment, and elevations of such 
interpreted features may differ by as much as .1 to .2 meters within a given pond. 
Compass bearings are probably accurate only to + 2-3 degrees, and taped distances may 
contain errors as large as +3% of taped distance due to interference by vegetation. The 
map positions of pond boundaries are based on a limited number of field measurements 
and in some instances are quite generalized. Overall, these data are sufficient for their 
intended screening and classification purposes but do not provide highly detailed 
topographic definition to the ponds. 


Data Processing and Pond Models 


Using the surface-plotting package SURFER (Golden Software, Inc., 1994), a 
kriging routine translated the elevation cross-sections and boundary positions into 
interpolated three dimensional models of the individual ponds. I then used the 3-D grids 
to calculate changes in wetted pond surface area (a measure of potential germination area) 
and pond volume (a measure of pond water balance) with change in the elevation of the 
modeled pond surface. The pond model output was then imported to the graphical 
plotting package GRAPHER (Golden Software, Inc., 1993) and displayed as stage/area 
and stage/volume response curves characteristic of the modeled geometry of the 
individual pond. The stage/surface area curves are a graphical representation of surface 
exposure with change in water level and can provide the basis for categorizing 
morphometric pond types. Stage/volume curves (not presented here) represent water 
storage within the modeled pond and may be used in conjunction with measured 
catchment areas to examine questions of pond water balance. 


Discussion 


Visual inspection of the interpolated pond grids shows considerable variation in 
basin shape and complexity. Figure 1 shows a selection of simple, nearly 
equidimensional pond models and their stage/area response curves. Ponds in this figure 
are shown at different scales; hypsometric response curves are at a common scale and 
show the contrast in pond size as maximum values along the area (horizontal) axes. 
Ponds of this general map pattern (ratio of major to minor axes < 1.5) are classified as 
‘equi’ in Table 1. 


Figure 2 shows a selection of elongate pond models with length to width ratios 
greater than 2, again shown with stage/area response curves. More complex pond model 
forms represent ponds with multiple basins, moat - like ponds with a central elevation 
maximum, ponds with lobate outlines, and one pond with an upland island within the 
wetland perimeter (Figure 3). In the latter set, Table 1 distinguishes between ‘irregular’ 
ponds (having complex outline but a single principal basin) from ‘complex’ ponds with 
multiple basins forming separate ponds under low - water conditions. 


The pond models reflect interpolated depths based on SURFER’s kriging routine, 
so measured depths at cross-section data points are not preserved exactly in the 3-D 
mesh. Basin shape and depth along measured cross-sections is reasonably well 
preserved, however. Artifacts of the kriging process are visible in many of the pond 
models, especially where there is uncertainty in the interpretation of the maximum pond 
water level. Upland areas above the interpreted full pool elevation are given artificial 
elevations of 0 meters to clarify the plotting of the model basins, therefore no reality 
should be ascribed to the “plateaus” surrounding the model basins. 


Figure 4 shows a composite plot of the depth - area curves for all 31 modeled 
ponds. This figure summarizes several basic pieces of information resulting from the 
pond surveys. Maximum model depths at full pool range from .5 m to 2.2 m, witha 
modal value near 1 m. Area of model “substrate” flooded at full pool ranges from less 
than 1000 m’ to almost 13000 m?. Overall steepness of the modeled pond bottom is 
represented by the average slope of the stage/area response curves. There is a visually 
apparent correlation between maximum flooded area and gradient of the substrate. Larger 
models (and ponds) have generally low gradients relative to smaller ponds. The two 
model ponds with the largest amplitude (EO 062 and EO 083) represent sites occupying 
relatively deep basins with topographic spill points well above the usual range in water 
levels. 


Statistical analyses of these curves are beyond the scope of this data collection 
project, but inspection of Figure 4 suggests that there are four geometric pond types in the 
surveyed population which can be defined by their hypsometric response functions. These 
four types should be represented in any set of sites chosen for analysis of the influence of 
groundwater or of the ecological effects of water balance change. Figures 5 through 8 
show these modes of model geometry plotted with their corresponding element 


occurrence numbers. Small, relatively steep geometries form the largest apparent cluster, 
identified as hypsometric class 1 in Table 1. Large, low gradient sites constitute another 
apparent mode (hypsometric class 3). Three relatively deep models of intermediate size 
and slope also stand out, and are designated class 2 in Table 1. Finally, the deep, high- 
amplitude basins discussed above form a distinct grouping, referred to in Table 1 as 
hypsometric class 4. 


Arduously collected field observations (the ‘wader penetration factor’) indicate 
that the large, low-gradient ponds have greater accumulations of organic sediment than 
most of the smaller ponds. In some cases more than a meter of peat underlies the central 
region of larger ponds, typically with H. aquatilis occurring along pond margins in 
slightly deeper, less infilled pools or ‘moats’. Survey elevations (and pond models) 
reflect this peaty surface. Pond contours and hydrologic response of these ponds would 
probably change with negative water balance changes, as organic sediments experienced 
greater surface oxidation. 


Table 1 shows specific conductance data for all sites visited twice during the 1997 
field season; Table 2 summarizes geochemical measurements for all ponds inventoried 
during 1997. Consistent with the results of Lesica (1992), specific conductance of the 
ponds varied from < 30 uS/cm to 400 pS/cm, with most ponds below 150 S/cm (Figure 
9). Measured pH ranged from 6.2 to 7.8, with most measurements between 6.5 and 7.5. 
pH is diurnally variable in dilute surface waters supporting photosynthetic activity. Since 
collecting measurements at a consistent time of day was not possible for this survey, little 
can be said regarding pond pH beyond the fact that the ponds do not depart very far from 
neutrality and that there is a general correlation between values of pH and specific 
conductance. Of the 34 ponds with early and late season specific conductance 
measurements in 1997, only five showed apparently significant increases in 
concentration, despite consistent reductions in pond volumes. Several ponds showed 
substantial declines in specific conductance with declining volume, while most remained 
unchanged within the probable precision of our field measurements. Above - average 
precipitation during the summer of 1997, plant uptake of solutes, groundwater outflow, or 
perhaps (in the most concentrated ponds) carbonate equilibria are possible explanations 
for this general lack of evaporative concentration with declining pond volume. 


Most H. aquatilis ponds have been considered topographically closed under 
present climatic conditions (Shapley and Lesica, 1997). Careful observation of the 34 
surveyed ponds showed that at least 12 have spill points occupied frequently enough to 
maintain some channel morphology. Outlet elevation control of pond hydroperiod and 
interpond exchange of surface water during wet periods thus appears to be more common 
than previously supposed. 


This project provided a broad sampling of pond morphomotry and defined 
variability significant to the expected interaction of H. aquatilis and hydrologic units of 
the species Montana habitat. Further analyses of pond morphometry should relate pond 
basin data to catchment size and topographic relief, likely to be important controls on 
shallow groundwater gradients and water balance dynamics. Water chemistry is 
generally dilute but displays a range of solute concentration that may reflect differential 
ground-water relations, therefore site selection for detailed hydrologic analysis should 
include measures of water chemistry. The presence of pond outflow channels is a 
measure of basin water balance and water level stability and should be retained as a 
screening criterion for study site selection. The role of pond size and orientation on light 
regime should be evaluated with respect to management of adjoining lands. 
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Figure 5. Hypsometric response class 1. 
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Figure 6. Hypsometric response class 2. 
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Figure 7. Hypsometric response class 3. 
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Figure 8. Hypsometric response class 4. 
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Figure 9. Specific conductance and pH of pond waters, July and August 1997. 
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29 | 7.7| 296 | 7 19 |present w/ submerg. fl. CIRARV @ N side 

30 | 6.8] 56 | 718]; present w/ emerg. fl. 

31 | 7.4| 225 | 7 19 |none found PHAARU @ S side 

32 |6.9| 119! 88 {30 w/emerg. fl+fr. on edge of open water _|TYPLAT w/ HOWAQU 


33 | 6.8| 143 | 88 |few plants w/ fr. in open water moat @ S tip 

34 L 6.8; 77 | 88 |few plants w/ emerg. fr. among sparse CARVES 
35 | 6.1} 38 | 88 |50 plants w/ emerg. fl. + fr. —__ algae bloom 
36 | 6.6 | 139 | 82 |20 plants scattered in open water @ SW side 

37 | 6.6| 145 | 83 |small plants w/ fr. among emergents @ S end ca 
3816.5) 62 | 83 (2 plants w/ emerg. fl. + fr. in open water nr S shore 


Table 2. Field pH and specific conductance of the inventoried ponds. 


[Exotics 


Howellia 


58 | 83 |20 plants w/ emerg. fl. + fr. in open water @ N side | 
none found; most of pond surveyed 


43 | 6.6| 65 [8 2 |none found; deepest not surveyed 
44/65) 58 | 82 [3 plants w/ submerg. fr. on edge of open water 
45 | 7.3| 93 | 84 |50+ w/emerg. fl. @S shore + frag. PHAARU in middle 
46 |6.9| 69 | 84 jnone found; survey confined to edges 

113 | 84 |none found 

126 | 84 |20 w/ emerg. fl. in deep open water 
49 | 7 | 98 | 88 (50 w/emerg. fl. + fr. in open water @w + S ends 


51/62] 31 | 84 [20 w/emerg. fl. in open water center in Sium 

52 | 6.3| 34 | 8 11 |<20 plants w/ emerg. + subm. fr. in Typha @ N end 
79 | 722 |few w/submerg. fr.@ W side + N end PHAARU @ N end+E side 
195 | 731 |50+ w/ fr. in open water + in CARVES 
152 | 731 |fragments along S shore 


56 | 6.6} 40 | 83 |<10 plants w/emerg. fr. @ S end in Sium 


57 | 7.3 | 288 | 728] none found on edges; deep center not surveyed 


58 | na Beeches See ar teeta 


59 | 6.7} 50 | 723 |none found; deep water not surveyed 

60 | 7 
61 | 6.8| 45 | 723 |20 w/ emerg. fl. in open water with Sium 

62 | 6.9| 44 | 723 jnone found; deep water not surveyed 

7 30 {none found; all but deepest searched 

few plants w/ emerg. fl. + fr. in open deep water 
65 | 6.9; 30 | 730 |none found along S shore 

66| 7 | 45 | 730 {few plants w/ fr. among heavy POTGRA @ S end 
67 | 6.1] 42 | 730 |fragments w/ fr. In open water @ N + S ends 

68 | 6.8! 35 | 730/50 w/emerg. fl.+fr. @ S side in open water 

69 | 6.9] 49 | 730 jnone found along margins 

70 |}7.5; 74 | 718 |few submerg. 

71 | 6.6| 110 | 724 |<10 w/submerg. fr. + emerg. fl. @ N end 
72 | 7 | 230] 84 |present w/ emerg. fl. PHAARU 
73 | 6.7| 124 | 724 /20+ on shaded edge + with Equisetum; center unsurveyed 

74 | 6.6) 94 | 724 |few w/emerg. fl. in shade @ NE side 
75 | 6.9! 180 | 724 |< 20 w/ emerg. fl. across N 2/3 of pond 
51 | 724 jabundant in Carex on N side 


—{__ 


PHAARU on E side 


Table 2. Field pH and specific conductance of the inventoried ponds. 


Howellia 


Exotics 


—{ 


81 |20+ w/ fr. in sparse CARVES + shaded W end 


8 1 |fragments in pond center 


7 25 | fragments w/ fr. in bay on S end 


725 none found; center unsurveyed 


7 24 |20+ w/ emerg. fl. in open water of sidearm 


210 


7 24 |20+ in open water; incompletely surveyed 


87 | 724 


none found 


;| 400 


7 25 |1 plant w/ emerg. fl. 


PHAARU; hvy algae bloo 


85 | 6.7| 65 


731 |few frag. w/ fr. 


7 | 194 | 


731 |20 w/ emerg. fl. @ E side; frag. @ W side 


7 | 220 


731 [50 w/emerg. fl. @ E side in sparse CARVES, EQUFLU 


137 | 


7 31 (none found; deepest unsurveyed 


731 |100 frag. w/ emerg. fl. + fr.; none rooted 


mod. algae bloom 


7 22 |none found; deepest unsurveyed 


7 22 |few w/ submerg. fl. @ open water in Sium 


[7.2] 72 


—t- 


7 22 |none found; deepest unsurveyed 


6.6; 89 
150 
173 
184 
111 
108 
7.1} 93 


66 | 723 


<5S0 w/ emerg. + submerg. fl. 
none found; deepest unsurveyed 

none found in complete survey 

none found; most of pond surveyed 

none found in complete survey 

several w/ submerg. fl. among logs @ S end 


abundant w/ emerg fl. + submerg. immat. fr. 
abundant w/ submerg. + emerg. fl. 
10 w/ submerg. + emerg. fl. @ N+S sides 


few w/ submerg. fl. @ W side open water in Sium 


CIRARV nearby 


<50 w/ emerg. + submerg. fl..@ S end among logs [ 


100+ w/ emerg. fl. 


50+ w/ emerg. fl. + fr.; sparse @ E+W sides 


Table 2. Field pH and specific conductance of the inventoried ponds. 


